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1.  Introduction 


Experimental  knowledge  of  the  flow  field  generated  by 
rotating  turboimpellers  is  essential  for  the  research  and 
development  of  turbomachinery.  This  information  is  used  to 
refine  design  methods,  develop  new  flow  models  which  include 
secondary  flow  and  tip  clearance  effects,  and  especially  to 
verify  computer  programs  designed  to  calculate  flow  through 
rotating  blade  rows. 

Laser  velocimeters  have  been  used  successfully  in  recent 
years  to  measure  the  flow  inside  and  downstream  of  rotors 
(see  Ref.  1) .  Certain  disadvantages  have  become  apparent, 
however.  The  laser  techniques  are  reliable  only  in  the  hands 
of  experienced  investigators,  the  pressure  field  remains 
unknown,  and  usually  the  measurement  of  more  than  two  compo¬ 
nents  of  the  velocity  field  is  complicated  and  expensive. 
Furthermore,  it  is  difficult  to  perform  measurements  close  to 
walls.  Development  of  alternative  techniques  to  overcome 
these  deficiencies,  as  well  as  to  achieve  redundancy  in  mea¬ 
suring  the  flow  field,  are  reasonable  and  worthwhile  tasks. 

This  report  describes  a  particular  method  and  the  com¬ 
putational  support  necessary  to  measure  the  flow  field  behind 
an  impeller  in  the  stationary,  bladeless  gap. 


2 .  Description  of  Method 

The  following  method  requires  two  semiconductor  pressure 
probes  along  with  a  technique  for  synchronized  sampling  for 
determining  the  fluid  velocity  vector  downstream  of  a  rotor. 

The  two  probes  (see  Fig.  1)  are  positioned  inside  the 
machine  casing  so  they  will,  in  turn,  intercept  periodically 
the  same  part  of  the  flow  leaving  a  particular  passing  rotor 
passage.  Each  probe  reading  is  sampled  when  the  designated 
blade  passage  reaches  a  desired  position  relative  to  the 
probe.  Synchronization  is  achieved  through  a  suitable  method 
(Ref.  2,3). 

Four  quantities  are  needed  to  determine  the  velocity 
vector:  yaw  angle,  pitch  angle,  static  pressure  and  total 
pressure.  Accordingly,  four  measurements  must  be  made  to 
evaluate  these  unknowns.  By  rotating  the  probes  about  their 
tips,  pressure  readings  in  four  different  directions  can  be 
taken,  and  the  data  used  to  calculate  the  velocity  vector. 
Computer  program  VELOCITY,  given  in  Appendix  II,  was  developed 
to  perform  the  somewhat  arduous  calculations. 

The  geometries  of  the  two  probes  are  shown  in  Fig.  1. 
Before  being  used,  the  probes  must  be  calibrated  so  their 
responses  to  flows  coming  from  different  directions  are  known. 
A  highly  directional  probe  is  desired  to  increase  the  accuracy 
in  finding  the  yaw  and  pitch  angles,  and  consequently  the 
velocity  magnitude.  The  following  method  is  recommended  for 
calibrating  each  probe  - 
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1.  Establish  a  steady,  controlled  flow  of  fluid,  and 
determine  the  velocity  vector  at  a  certain  region  of 
the  flow. 

2.  Position  a  probe  in  the  flow  and  rotate  the  tip  so 
that  a  sequence  of  pressure  readings  are  taken  for 
a  constant  yaw  angle  and  a  varying  pitch  angle. 
Repeat  the  procedure  at  a  new  yaw  angle  using  the 
same  pitch  angles.  The  result  will  be  an  array  of 
pressure  readings  corresponding  to  a  set  grid  of 
yaw  and  pitch  angles  (Fig.  2)  . 

3.  From  the  known  flow  velocity  and  pressure  readings, 
a  coefficient  of  pressure  can  be  calculated  for  each 
angle  set: 

p-ps  C  =  Coefficient  of  pressure 

C  =  — —  where :  p 

p  PT  ps  p  =  pressure  reading 

pg  =  static  pressure  of  flow 

PT  =  total  pressure  of  flow 

The  table  of  C  's  as  well  as  the  yaw  and  pitch  angles  which 
P 

correspond  to  them  are  now  in  the  form  required  for  input  to 
program  VELOCITY. 

The  probe  calibrations  should  be  insensitive  to  Mach 
number  and  pressure,  and  are  not  valid  for  supersonic  flows. 
Should  any  significant  variations  in  C^  be  observed  for 
different  flow  conditions,  further  calibrations  will  be  re¬ 
quired  and  an  additional  iteration  scheme  added  to  the  com¬ 
puter  program. 
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Experience  with  the  two-probe  technique  has  shown  that 
excellent  results  are  achieved  when  a  probe  type  A  is  rotated 
to  the  three  positions  +25°,  0°,  -25°  yaw  at  0°  pitch) (and 
probe  type  B  is  used  at  0°  yaw  and  25°  pitch.  Fig.  3). 

The  two-probe  technique  is  strictly  applicable  only  to 
periodic  flows.  However,  data  obtained  on  successive  rotations 
of  the  rotor  can  be  averaged  to  eliminate  non-periodic  fluctua¬ 
tions.  This  was  effective  for  tests  reported  in  Ref.  2.,  where 
a  single  probe  was  used  to  establish  the  peripheral  blade-to- 
blade  distribution  of  flow  yaw  angle. 

It  is  noted  that  the  method  reported  here  is  a  further 
development  of  that  reported  earlier  in  Ref.  6,  and  overcomes 
some  of  the  earlier  limitations. 
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Figure  3.  Orientation  angles  of  the  probes  relative 
to  the  laboratory 
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3.  Theory 


The  velocity  vector  for  a  three-dimensional  flow  can  be 
described  with  three  scalar  quantities.  The  nature  of  the 
problem  suggests  using  two  angles  (a  yaw  angle  and  a  pitch 
angle) ,  and  the  magnitude  of  the  velocity  (Fig.  4) . 

Since  pressures  and  not  the  velocity  are  measured,  the 
static  and  total  pressures  must  first  be  determined,  and 
Eq.  (1)  used  to  evaluate  the  velocity. 


£t 

?s 


(l  +  Iji  mV/y-1 


(l) 


Altogether,  four  unknowns  need  to  be  evaluated:  the  yaw  and 
pitch  angles,  and  the  total  and  static  pressures. 

Four  equations  are  needed  to  determine  the  four  unknowns 
They  are  derived  from  the  four  pressure  readings,  each  pres¬ 
sure  reading  having  been  taken  in  a  different  direction  as 
described  above.  The  following  equations  for  the  coefficient 
of  pressure  can  be  written: 


1  ps 


PT'PS 


i  =  1..4 


(2) 


The  C  .'s  are  a  function  of  the  orientation  of  the  probe 

pi 

relative  to  the  flow;  i.e.  ,  for  a  given  flow  the  measured  C 


will  vary  measureably  as  the  probe  is  turned  into  and  away 


<£-  PITCH  AiiGLc 

V-  IIVII-  MAGNITUDE  OF 

VELOCITY  VECTOR 

X -  REFERENCE  FRAME  FIXED 
IN  THE  LABORATORY 

•  ► 

Figure  4.  Velocity  Vector  V 


from  the  flow.  Each  "probe"*  will  have  its  own  C  charac- 

P 

teristics  determined  experimentally.  The  result  will  be  a 
table  of  Cp  vs.  yaw  and  pitch  angles  for  each  probe. 

cpi  =  function  (aRi  ,  $Ri)  i  =  1..4  (3) 

For  realistic  problems,  only  one  point  (a  ,  <p)  exists 
where  the  C^'s  in  Eq.  (2)  will  equal  the  of  Eq.  (3)  for 

the  four  probes'  pressure  readings. 

The  probes'  characteristics  (C^'s)  are  in  tabular  form 
because  they  cannot  be  represented  analytically  due  to  the  stem 
effect  and  production  inaccuracies.  Therefore,  a  numerical 
solution  to  the  problem  is  required.  The  procedure  chosen  for 
solving  the  problem  is  a  systematic  trial-and-error  search 
process,  essentially  a  convergence  scheme  on  two  variables: 
yaw  angle  and  pitch  angle. 

The  flow  direction  is  assumed  to  fall  within  some  set  of 
bounds,  defining  the  search  area  for  yaw  and  pitch  (Fig.  5). 

By  setting  up  a  grid  of  points  in  this  region  and  checking  how 
well  each  point  satisfies  the  criteria  of  equality  of  coeffi¬ 
cients  of  pressure  (C  Vs)  calculated  with  Eqs.  (2)  and  (3), 
the  point  with  the  smallest  error  can  be  found  and  used  as  a 
first  approximation  to  the  solution.  Repeating  this  procedure, 
only  with  a  smaller  grid  and  search  region,  will  result  in  a 
better  approximation.  This  sequence,  represented  in  Figs.  6 


*  Here,  the  term  "probe"  refers  to  a  particular  probe  type  in 
a  particular  position. 
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PITCH 


*  FLOW  DIRECTION  OF  THE  FLUID 
-  BOUNDARY  OF  SEARCH  AREA 


Figure  5.  Search  area 
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and  7  is  repeated  until  either  the  desired  accuracy  is  reached 
or  fatigue  sets  in.  Program  VELOCITY,  described  in  the  follow- 
mg  section,  was  written  to  perforin  these  calculations. 
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Figure  7.  Flow  Chart  of  the  Search  Procedure 


4.  Program  VELOCITY 


Program  VELOCITY  was  written  to  perform  the  calculations 
outlined  in  the  previous  section.  A  description  of  the  program 
and  its  subroutines  is  given  below.  Fig.  8  surnarizes  the 
major  sections  and  organization  of  the  program. 

For  each  run,  program  VELOCITY  reads  the  calibration 
tables  for  the  two  probes  from  files  outside  the  proqram. 

(Input  formatting  is  discussed  in  Appendix  V.)  Subroutine 
INPUT  performs  the  necessary  work,  and  can  be  modified  to 
accommodate  different  input  schemes  if  desired. 

The  fluid  temperature  and  molecular  weight  are  entered 
next.  These  properties  are  assumed  to  remain  constant  through¬ 
out  the  run. 

The  settings  for  each  pressure  reading  are  read  next.  A 
setting  contains  the  following  data:  probe  type  (A  or  B) ,  yaw 
angle  setting,  and  pitch  angle  setting.  Again,  these  settings 
will  not  change  for  the  duration  of  the  run. 

Finally,  the  four  pressure  readings  are  entered. 

The  first  scan  is  initiated  and  covers  the  entire  region 
of  expected  flow  directions,  -40°  to  +  40°  in  both  yaw  and 
pitch  angles  in  the  present  case.  Points  are  chosen  every  5°, 
each  point  representing  a  unique  pair  of  yaw  and  pitch  angles. 
For  each  point,  a  static  pressure,  a  dynamic  pressure,  and  an 
error  are  calculated  by  the  scheme  described  below. 

A  point,  say  (a  ,  <M  is  tested;  i.e.,  a  test  is  performed 
to  prove  whether  assumed  flow,  oriented  a  degrees  yaw  and  * 


degrees  pitch  relative  to  the  laboratory  reference  frame, 
corresponds  to  the  four  pressure  readings.  The  direction  of 
the  flow  relative  to  each  probe  setting  is  calculated.  For 
probe  setting  i  ,  oriented  at  (ct^  ,  relative  to  the 

laboratory,  the  assumed  flow  approaches  at  a  relative  angle 
of : 


4> 


Ri  = 


(5) 


where  (otRi  >  ‘t1^)  are  the  Vaw  anc^  pitch  angles  respectively 
of  the  assumed  flow  relative  to  probe  setting  i  .  The 
calibration  table  for  the  probe  used  in  setting  i  is  con¬ 
sulted  and  a  C  (a  .  ,  <)>  . )  returned.  Subroutine  CPCAL  locates 
p  K1  K1 

or  calculates  the  desired  C  values  in  the  table.  The  scheme 

P 

used  in  CPCAL  is  a  search  technique  to  find  the  values  of  yaw 
and  pitch  surrounding  the  desired  point,  and  then  a  linear 
interpolation  over  these  four  points  as  shown  in  Fig.  9. 

Eq.  (2)  can  be  rewritten  in  the  form 

(CPi)PT  +  (1'Cpi)PS  =  Pi  i=1*-4  (6) 

the  only  unknowns  being  pT  and  pg  .  With  four  equations  and 
two  unknowns,  the  problem  will  be  inconsistent  unless  the  true 
a  and  <t>  were  chosen.  Accordingly,  the  following  schemes 
were  used  to  evaluate  p  ,  p„  and  an  error. 

S  1 


17 


Firurp  Q.  Li'p,--  i  o  terpol  at  ior  betv;r?r,  four  noir.ts  to  finl 


Define : 


C  =  £  C 
-£  i=i  pi 


E  ■  Si 


C  =  minimum  of  (C  ,  C  ,  C  ,  C  ) 

pm  Pi  p2  P3  P4 


Pm  =  corresponding  to  the  C  chosen  above. 


(Cp)pT  +  <4-Cp)ps  =  P 


and  also 


<C  >PT  +  (l-c )ps  -  P„ 
rm  m 


These  two  equations  can  be  solved  for  pT  and  p 


P(l-C  )  -  p  ( 4-C  ) 

*  Pm  m  -£ 


PT  -  c  -  4C 
P  Pr 


c  (P  )  -  C  (p) 

-E  m  PmC 


C  -  4C 
P  Pr 


‘i  p  .  ~p 

Error  =  L  ABS(C  -  -  ~  -s-)  /4 
i=l  pi  Pt'Ps 


These  schemes  were  chosen  for  two  reasons: 

1)  They  used  all  the  available  data  to  derive  an  error 
which  would  effectively  represent  the  accuracy  of 
the  guess. 

2)  No  singularities  in  the  calculations  can  occur  except 
for  the  case  of  four  equal  C^'s  (which  physically 
represents  trying  to  find  an  intersection  point 
among  four  parallel  lines).  If  the  measurements  are 
taken  in  the  suggested  directions,  this  anomalous 
point  will  not  appear. 

For  each  point  guessed  in  the  initial  scan,  an  error  is 
calculated  and  the  point  with  the  smallest  error  is  saved.  A 
new,  finer  search  grid  is  composed  using  this  point  as  the  new 
origin.  The  boundaries  of  the  new  grid  are  the  points  from 
the  old  grid  which  were  closest  to  this  new  origin.  Referring 
to  Figure  10,  if  x  represents  the  true  solution,  the  new 
boundary  would  be  formed  by  the  points  marked  B-I,  and  the  new 
grid-width  would  be  one  third  as  large.  This  factor  was  chosen 
to  minimize  the  number  of  guess  evaluations.  (The  first  scan 
contains  a  large  number  of  guesses  in  order  to  correctly  isolate 
the  general  region  of  the  solution) . 

The  search  procedure  is  performed  on  each  new  grid,  and 
the  process  repeated  until  the  grid  width  is  less  that  0.5°. 
After  the  final  scan,  the  best  guess  is  used  to  calculate  the 
flow  velocity  and  Mach  number.  The  results  are  printed  out 
and  the  next  four  pressures  requested.  If  no  values  are 
entered  (end  of  data  set) ,  the  program  ends. 

20 
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•*  TRUE  SOLUTiOr: 

fi.jure  10.  Defining  new  grid  b&ur.duriei  fron  tne  nearest  neighbors 
of  the  point  with  the  smallest  error 
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5. 


Discussion 


Extensive  tests  with  program  VELOCITY  have  led  to  the 
observations  and  suggestions  listed  below: 

1.  Excellent  results  are  achieved  when  the  probe  settings 
are  at  (yaw,  pitch)  angles  of  (-25,0),  (0,0),  (25,0) 
and  (0,25)  degrees.  This  corresponds  to  a  rotation 

of  probe  type  A  from  -25°  to  0°  to  25°,  and  one  read¬ 
ing  from  probe  type  B  at  (0,25).  Poor  results  were 
achieved  for  the  symmetric  case  of  readings  at  (+25,0) 
and  0,+25)  degrees. 

2.  Highly  directional  probes  increase  the  accuracy  of 
the  procedure,  especially  if  the  variation  is 
significant  when  the  flow  is  nearly  head-on  To 
achieve  these  characteristics,  the  following  design 
suggestion  is  offered.  The  probe  can  be  formed  with 
a  spherical  tip,  the  pressure  tap  being  located  in 
the  center.  To  prevent  damage  to  the  sensitive  trans¬ 
ducer  located  behind  the  pressure  tap  and  to  improve 
the  frequency  response,  the  void  between  the  pressure 
tap  face  and  the  transducer  should  be  filled  with  an 
an  appropriate  liquid  and  the  opening  of  the  pressure 
tap  sealed  with  a  thin,  low-inertia  membrane. 

3.  Higher  accuracy  naturally  results  if  more  calibration 
points  are  taken  for  the  probes'  tables.  The  linear 
interpolation  scheme  can  be  replaced  by  the  second 
order  scheme  offered  in  Appendix  5  (if  no  significant 


22 


anomalies  occur  in  the  calibrations) ,  the  second  order 
method  requiring  fewer  calibration  points  (say  every 
15°)  than  the  linear  method  (every  5°  or  10°). 

The  use  of  two  probes  of  relatively  simple  geometry 
in  periodic  flow  is  less  cumbersome  and  complex  than 
the  use  of  five-hole  probes  (Ref.  4). 
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Notation  Summary 


c 

p 

Coefficient  of  pressure 

C  is  a  function  of  a  and  <!>  ,  C  =  C 

P  P  P 

cpi 

- 

Coefficient  of  pressure  for  probe  setting 
Cpi  =  Cp(aRi  '  ^Ri^ 

- 

Siam  of  the  four  C  .  's 

P1 

cpm 

- 

Minimum  of  the  four  C  .'s 

Pi 

Cpn 

- 

cp<“i  1  *1* 

Cp12 

- 

Cp(al  '  *2> 

C 

P21 

- 

Cp(0,2  >  *1’ 

C 

p22 

- 

cp(a2  -  *2> 

P 

- 

Pressure  (all  pressures  are  absolute) 

P. 

1 

- 

Pressure  read  from  probe  setting  i 

ps 

- 

Static  pressure 

PT 

- 

Total  pressure  (stagnation  pressure) 

E 

- 

Sum  of  the  four  pressures  (Pi's) 

m 


Pressure  at  the  setting  where  C  occurred  (i.e., 

pm 

P  =  P.  ,  where  i  =  m  ,  defined  in  C_  ) 

m  1  Pm 
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V  -  Velocity  magnitude  of  the  fluid  particle 

V  -  Fluid  velocity  vector 

a  ,  <J>  -  Yaw,  Pitch  angles 

(*1  ,  <|>i  -  Yaw,  Pitch  angles  for  probe  setting  i 

aRi  ’  ^Ri  "  Yaw'  angles  for  the  assumed  flow  direction 

direction  relative  to  the  probe  setting 

p  -  fluid  density 
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'opendix  I  -  Program  VELOCITY 


VPLCUTY  A  NO  DIRECTION 
Gl  A  'LUIL  AT  A  POINT  UGlHC  f  rur  pressures 


~A.J,_  TAYl,_K 


.ivy  thf  c  fuocs  srs.iLD  ev  pr.cc? so  in  feu* 
distinct  niFrairK:,  a  no  kwain-;  th-;  characteristics  cf  if 

I’HClIt  I  AN'  FLU  10  PRCPCRT  ICS.  THC  FLLD  VEICCITy,  riPI'f  Tfi'N. 
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.-ALL  INPIfNALPF  l  •  NPH  I  1  ,  RRH  I  .CP1.IEKI 
CLL  1  NP  T  (N  Al.  PE?  .  NPI11  2  ,  PPiJ2  .CF2  .  I  CP  > 

I  r  (  i ;  . .  nf  .  o »  ;t  t-.:  i  c?o 


-ia.  n-r  M.LtFMNu  ri.uic  pfchftics: 


»v 

GA  K  A 
TC. 

(.U  N  > 


YCLLCULAP  i'FjHI  cr  THE  FLU  10 
14TIC  cr  FPtCIFIC  HEATS  CF  THC  FLU  T  0 

temper  atlm:  elgrels  ccnticraoc  cr  iff  c  a  3 

*:  i  T  I  A  A  1  r  rf  TfE  Cl,  Mp|<f.  FSICIUITV  FACTOR 


50  ,i  £  A  Z  (  ■:  .  5  C  I  0  )  «  N  ,  0  AYVA.TC.CCVC 

5  0  10  F  'J  fi V  A  J ( 4  F  1 C  .  4 ) 

rgaf  -  f:i4./v»n 

*»IK(?|77CC)  ¥M  ,C.  A  v,4A  .  TC  .  CLMF 

7/ 00  FCliKAIC  FL  IC  FOCPcHTICS  :*.//.*  MOLECULAR  »T  =  •  .  T  ?0  .F  3 . 4  ,  /  , 

1  •  FATIC  CF  OPCCIFIC  FEATS  =  •  ,T30  .F<?.4  ./ .  •  TEMPERATURE*  * 

;  *  :cc  C  =  *  .  T JO ,T c .4 . / ,  •  CCVP  (CSSIEILTY  FACTrfi  =  *  .  T3C  .F8.4 ) 

am  IE  (  C  .  t  CC  C  ) 

G)C0  -:C.H'4I!M  STATIC  T  CT  AL  *  .  I  2K  .  •  V  A  *  PITCH*.:*, 

1  'VELOCITY  MACH *,/,  3X ,* PRESS  (PA)  PPCSS(PA)*, 

3  P* , 'ANCLE  ANGLE • .6X, • ( M/SCC I  NUM2CS*.//) 

«  P  I  T  :  <7 , 7  7  1  C  ) 

7710  FUR  NA  T  ( // / ,  •  PPr-lE  YAW  PITCH  FREFiURE',7. 

I  •  I Y  PE  SETTING  SETTING  PC  A3  (PA  )*,//> 


«*«**  FT  A  1  ,*  T  L  r  C  P  («*♦♦ 


Hi_  A,; 


in  the  experimental  jata  f  fjp  this  ot  tt  pv  in  a  t  ic  n 


np  r-e  ( i  ) 

ALM  I  ) 

PH  I  <  1  ) 

ppe:  s(  I ) 


TFT  PPCI’C  TYPE  (EITHER  t  CR  2)  CF  PPCLF  SETTING  I 
ALPHA  (YAW)  ANGLE  QF  PROBE  JETTING  I 
PHI  (PITCH)  ANGLE  EF  PPEBF  SETTINC  I 
PF'-FGUFE  REAP  l)Y  PRCOE  SETTING  I 


ALM.  c  H I  I  ,  NPNOI  .  CONTAIN  THE  NEW  VALUES  CF  THE  YAW,  PITCH. 

ANC  PxC’P  TYTE  TCP  EACH  SETTING.  IT  THC  VALLc  READ  TCP  THC  PPC 3E 
TYf_  1?  JEM.  IMV[II=0),  Tr-LN  THC  PHCPC  SETTINC  TCP  THE  °REVICLS 
IF  li  USE.:.  r.C  ’DEFAULT  VALUES  are  PRCVIGFC.  SC  THE  FIRST 

tr  i a i  *m:t  ccntain  the  probe  settings. 


10  ,  "  I  =  1  ,  4 

Jl  4  F  (  **  .  SC  J  IJ  |?NC=  MS  )  PRESS  (  I  )  .  ALP  I  .PHI  I  ,N  M  !  I 

10  30  -  l  o  .«  A  T  (  3  r  I  p  .  4  .  I  1  ) 

IMNT'HI  .'•;.))  CC  TC  15 
A  L  P  (  I  )  =  A  L  p  I 
•>F  I  (  I  >=PF  I  I 
,r,;i  t  i  >  =  n r i.  >  i 

15  A. II  (7.773?)  N.M'M  I )  ,ALP(  I  )  .PHI (  I  ) ,PHCS3(  l  ) 

7  7  3  0  ** C  r  v  a  t  (  i  >  ,  t  -  ,  7  >  ,  ?r  I  o  .  3  .  F  I  4  .  P  ) 

3  0  SCN  I  I N^E 

A  PI  I  L  ( 7  . 7  7  J C ) 

7  7  JO  •  Y  A  T  <  /  /  ) 

FSMIlIJF  rC  'IMNf  RANGE.  C1.IO  WICTH.  ANO  INITIALISE  ESPCR 
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0  0003  0  1  0 
C  C  0  C  0  0  ?  0 
CC000030 
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:  o  )  o  o  o  s  o 
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VELOCITY 
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!0t  L  - 

= 

A  V  «  »  =  4  C  . 
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F  V  A  >  =  ft)  . 

2  2  f  =  5  • 
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c 

V 
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TCh  ANCIC 

CL  F 
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c 
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APPENDIX  II 


VELOCITY  NOTATION  SUMMARY  -  main  program 


ALP (I)  -  Yaw  angle  of  probe  setting  I 

AMIN,  AMAX  -  define  the  minimum  and  maximum  yaw  (alpha)  angles 
of  the  search  grid 

COMP  -  the  compressibility  factor  of  the  fluid 

CP (K)  -  C  interpolated  from  the  appropriate  calibration  table 
for  p  probe  setting  K 

CPMIN  -  stores  the  minimum  C^  found  during  this  guess 

CPSUM  -  stores  the  sum  of  the  four  C  's  read  by  Subroutine  CP CAL 

P 

CPl ,  CP2,(I,J)  -  C  calibration  table  for  probes  1  and 

P 

C0  -  sonic  velocity 

DEL  -  search  grid  spacing  (degrees  of  angle) 

DENOM  -  stores  an  intermediary  mathematical  quantity 

ERRMIN  -  stores  the  minimum  error  found  so  far  for  the  problem 

ERRR  -  C  average  error  characteristic  for  the  guess 
P 

GAMMA  -  ratio  of  specific  heats  for  the  fluid 

IER  -  input  error  flag  =  0  means  no  error,  =  1  an  error  occurred 

while  reading  in  the  C  calibrations 

P 

IFL  -  interpolation  error  flag  =  0  interpolation  accomplished 

=  1  range  of  the  calibration  table 
was  insufficient 

MACH  -  fluid  Mach  number 

NALPH1 ,  NALPH2  -  number  of  yaw  angles  across  the  edge  of  the  Cp^, 
Cp2  calibration  tables 

NPHIl,  NPHI2  -  number  of  pitch  angles  across  the  edge  of  the  Cpl , 
Cp2  calibration  tables 

NPRB(I)  -  probe  type  for  probe  setting  I  (either  1  or  2) 

PHI  (I)  -  pitch  angle  of  probe  setting  I 
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PMIN ,  PMAX  -  define  the  minimum  and  maximum  pitch  (phi)  angles 
of  the  search  grid. 

PRB1,  PRB2  (N,J)  -  contains  the  alpha  and  phi  angles  for  use  with 
Cp^ ,  Cp2  respectively.  J  =  1  refers  to  yaw  angles 

J  =  2  refers  to  pitch  ancles 

PRESS (I)  -  pressure  read  at  setting  I 

PRMIN  -  stores  the  pressure  at  the  setting  corresDonding  to  CPMIN 
PRSSUM  -  stores  the  sum  of  the  four  input  pressures 
PSS  ~  contains  a  static  pressure  characteristic  for  this  guess 
PTT  -  contains  a  total  pressure  characteristic  for  this  guess 
RGAS  -  ideal  gas  constant  ( Joules/kg-°K) 

RHO  -  fluid  density 

TC  -  fluid  temperature  °C 

VEL  -  fluid  velocity 

WM  -  molecular  weight  of  the  fluid 

X , Y  -  yaw,  pitch  angle  guess  (one  of  the  search  grid  points) 

XMIN ,  YMIN  -  yaw,  pitch  angle  where  the  smallest  error  was  found 

XR,  YR  -  yaw,  pitch  angles  of  the  guess  relatwe  to  the  probe 
setting  being  considered. 


NOTATION  SUMMARY  -  SUBROUTINE  CPCAL 


AP, AN  -  Yaw  angles  above  and  below  the  desired  yaw  angle 
CP (NA , NP )  -  calibration  table 

Cll ,  C12,  C21,  C22  -  0^  values  surrounding  the  desired  C^ 

I FLAG  -  error  flag  =  0  means  the  interpolation  succeeded 

1  the  range  of  the  table  was  too  small 

MNA ,  MNP  -  Stores  the  location  of  the  calibration  yaw  (alpha) , 
pitch  (phi)  angles  below  the  desired  yaw  and  pitch  angles. 

MXA,  MXP  -  Stores  the  location  of  the  calibration  yaw,  pitch 
angles  above  the  desired  yaw  and  pitch  angles. 

NA,  NP  -  number  of  yaw,  pitch  angles  in  the  C  calibration  table 

P 

PP,  PN  -  Pitch  angles  above  and  below  the  desired  pitch  angle 
PRB (N,K)  -  Contains  the  yaw  and  pitch  angles  for  the  calibration 


table 

X, Y  -  Yaw  and  pitch  angles  where  a  C  is  sought 

XB,  YB  -  Fractional  distance  of  the  desired  yaw,  pitch  angle 
between  the  known  calibration  angles 


7  -  the  interpolated  C  value  for  X,  Y 
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NOTATION  SUMMARY  -  SUBROUTINE  INPT 


CP(I,J)  -  Calibration  table  read  from  the  file 


NA  -  Number  of  yaw  angles  on  the  edge  of  the 


C 

P 


table 


NP  -  Number  of  pitch  angles  on  the  edge  of  the 


C 

P 


table 


PRB (N , K )  -  contains  the  yaw  and  pitch  angles  for  the  C 
table  p 

K=1  yaw  angles 
K=2  pitch  angles 
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calibration 
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APPENDIX  m  -  Sample  Input 
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1  u  787  ij  , 

i  I!  VI  v  0  . 

.1  0  9 1  v  o  , 
1  0  S  9  8  0  , 
vfll  0  0 
t.  OS 4 SO  , 
1  0  8 '2(1, 
99)80 . 

).  2  o  a  o  o , 

i  I  S  0  9 1)  , 

1  1  S  0  1  (i  , 
.i. ;  '  0  j  (  i  , 

I.  )  884  0  . 

.1  ,1,  1)64  0  , 

1 1 668 0 . 

II.  1. 47  0  , 


i-40  £0,0  5,0 

~7S,0  0.  5 

o  •  0  . 

.  o  o .  ■■ 

11  ■  -:,S  .  (! 


0  •  0  .  i. 

0  ■  P.S .  0 

-as .o  o  i 

.?  S ,  0  o ,  -j 
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APPENDIX  ~J5L  '"Sample  Output* 


FLU  II)  FKiim  K  7  ILj  : 


1  U  t  u  «- 

UL  Ak  if,  T 

- 

: <  .hccc 

t  a  r  i  ; 

if  F  *  j)  i. 

c  1 1  i  c  r  l  a  i  :  - 

l  •  a  c  c  c 

rtf  v  v; 

h  4TU')r 

ICC  r  = 

:  l  .  )  c  c  o 

31.1  3  I  L 

T  y  K  ACT  L  - 

1  .  C  1  0 

DU 

>  a* 

J  I  I  Mi 

(••<1  1'J’JN;: 

IV  'I. 

it  [  >s . 

S  1:  1  1  1  y  o 

'■  (  a  (  , )  A  ) 

1 

-  2 5  .  C  J 

Vy  •  0 

It/tty 

.00 

1 

'1  .0 

c*  •  - 

1  C  C  V  c  c 

.0  0 

1 

2  1.0) 

l  .  r 

i  c.-c •: r 

*  ^ 
•  ‘  ~J 

C 

C  •  Vy 

~ .  co 

1  c  7  1  c  ( 

.00 

1 

: .  o  j 

o .  •: 

1  0  o  1  <  c .  c  0 

1 

C  .0 

c .  : 

loo  ? ye  .co 

1 

.‘z : 

c .  : 

! C70CC.CC 

j) 

C  .0 

2  : .  r  ) 

1  C  2  C  C  0 . 0  0 

1 

-Cl. 0 J 

i.  .  c 

1  c  ■>  3  J  c  .  c  0 

l 

C  .  3 

0  .  c 

!  C  •  •  t  c  C  .  C  0 

1 

2  e  . :  j 

: .  c 

lC0Ilt.ee 

i.' 

c  .o 

j  i .  c  0 

105 0CC  .30 

l 

,ao 

'  -i  1  c.'  .CO 

1 

C  .0 

A 

y  • 

1  C  5  5  0.00 

l 

o  i .  c 

C  .  C 

i  c  -j  i  r  .  3  c 

p 

c  •  0 

C  5  .  f  0 

‘  C  1  €  C  .  3  0 

1 

C  .  J 

c .  c 

l 

t  3  i>  C  y  •  )3 

j 

C  •  0 

:  .  .  C  9 

1 

l  f  •  0  <  C  .  0  c 

1 

J  J 

c .  c 

1 

r.:-  c  i  c .  c : 

1 

lL  -  *00 

0  .  0 

1 

1  2  C  1 C  .  C  0 

l 

0  •  j 

*J  •  c 

i  i  ?.  e  <  c .  c  o 

I 

C  .  0 

2  5.C0 

1  1  0  5  A  C  .  3  0 

l 

-21.00 

C  .  C 

1  IttEO.OC 

1 

21  .0  ) 

O  .  C 

1  1  1  3  7  J  .  0  C 

jIUU 
PftLij  (PA) 

re  r  t*. 

p  ;r_  js  c5  a  i 

V  A  Vt 
ANClt 

r1  l  T  £  F 

AN  CL  Z 

VlLClJ  TV 
( ►'/SEC.  ) 

»/  ACF 

M.WCCS 

i c c  ;o .'it 

i  1  C  1  1  £.  «  71- 

-  .  c .  7 

2  0.07 

l 2rt.7< 

C  .2  74 

1)0)5')  •  5  - 

1  10  C  s  <  .  1  0 

-•  •? .  i  c 

i . :  t 

1 1.  R  .  j  J 

C  .  1  72 

;o  j 19 .  i  > 

ll012V.lt; 

12.50 

0.3  7 

1H  7 .4  2 

C  .  1  A  4 

£-ifcCd.OJ 

l  lOUt.2- 

2  2.  15 

• lo  *ti 

1  C  )?  .  J  A 

C  .  1  a  5 

•i  t  e  e  <*  •  lo 

1  1 . 9e (  .  V  j 

-c.cc 

-0  •  J  U 

220.07 

c .  c .  t 

e  9L  79  .  o  -* 

120135.25 

-  7  .  C  .4 

.(.  2 

2  2  L  .  ? » 

C  .f  e  7 
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APPENDIX  V 


NOTES  ON  THE  USE  OF  VELOCITY 

INPUT:  The  required  input  consists  of  probe  calibration  data, 

fluid  properties,  and  finally  the  experimental  pressures.  Sub¬ 
routine  INPUT  reads  the  calibration  data  from  each  probe  type 
in  the  following  form: 

1.  The  first  card  contains  the  number  of  yaw  and  pitch  angles 
on  the  axes  of  the  calibration  table  (format  214) 

Ex:  19  19  means  19  yaw  and  19  pitch  angles  were  used  in 

the  calibration  and  the  C  table  will  therefore  be 

P 

19  x  19  in  size. 

2.  The  next  few  cards  contain  the  values  of  the  yaw  angles 

where  calibration  points  were  taken  in  the  table. 

Values  are  entered  in  format  F8.2,  one  angle  every  8  columns. 
After  all  the  yaw  angles  have  been  read,  the  pitch  angles 
are  entered  starting  on  a  new  card. 

3.  The  experimentally  determined  Cp's  °f  the  calibration 
surface  can  now  be  read  for  each  angle  pair  starting  from 
the  smallest  yaw  and  pitch  angle  and  with  the  pitch  angle 
varying  most  rapidly.  Ex. 

C  (-90),  -90),  C  (-90,  -80)...  C  's  are  read  format  F  8.5. 

P  P  P 

All  of  the  calibration  data  are  read  on  Machine  Unit  8: 

Cards  are  assumed  to  be  80  characters  in  length. 

The  following  fluid  properties  are  entered  next: 

Molecular  Weight 
Ratio  of  Specific  Heats 
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Fluid  Temperature  Deg  C 

Compressibility  Factor 

Machine  Unit  6  reads  this  data  from  one  card.  Format  4  F10A. 

At  last  the  experimental  results  are  entered.  Four  cards 
are  required  for  each  trial,  one  card  per  setting.  For  format: 

Columns  1-10:  Experimental  pressure 
11-20:  Yaw  Angle 

21-30:  Pitch  Angle 

31:  Probe  Type  (1,2,  or  blank) 

If  Column  31  is  left  blank,  only  the  experimentally  read 
pressure  is  registered;  yaw  and  pitch  angles  for  that  setting 
remain  unchanged  from  the  previous  trial.  The  first  trial  must 
contain  angle  settings  and  probe  type  since  no  default  values 
have  been  assumed.  Again  machine  unit  6  is  used  to  read  this 
data.  When  no  more  experimental  pressure  data  is  available, 
the  program  terminates. 

The  experimental  pressures  can  be  based  in  any  absolute 
system  of  measurement;  ex.:  Psia,  KPa,  Atm,  mmHg,  with  the 
same  numerical  results  (the  units  in  the  titles  of  the  static 
and  total  pressure  columns  will  not  apply) .  The  analysis  below 
shows  that  in  determining  velocity,  the  pressure  units  cancel. 
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The  velocity  is  calculated  from  V  =  MCq  ,  where,  from  Eq.  (1), 

1/2 


and  C  =  /yRT 

o 

Here , 

CQ  =  sonic  velocity 
M  =  Mach  number 
Ps'PT  =  static,  total  pressure 

R  =  ideal  gas  constant 


MW 

T  =  Fluid  Temperature  °K 
V  =  Fluid  Velocity  (m/sec) 
p  =  Fluid  density 
y  =  ratio  of  specific  heats 

CPCAL;  A  linear,  double-interpolation  scheme  is  employed  to 
determine  a  value  of  Cp  between  four  points.  A  second-order, 
double-interpolation  scheme  has  also  been  devised  and  tested, 
and  is  presented  at  the  end  of  this  report.  Figure  V-l  is  a 
graph  of  the  accuracy  of  both  schemes  as  a  function  of  the  num¬ 
ber  of  calibration  points  in  the  C  table.  Values  were 

It 

determined  by  filling  a  calibration  table,  extending  from 
-90°  to  +  90°  in  yaw  and  pitch  with  the  Cp's  which  would  result 
from  an  ideal  probe,  and  testing  6084  points  (78  x  78)  within 
the  table.  If  no  highly  unusual  distortions  in  the  calibrations 


of  the  probes  occurs.  Figure  V-l  shows  that  a  significant  reduc 
tion  in  the  amount  of  calibration  required  is  possible  with  a 
second  order  scheme.  Further,  if  the  accuracy  of  the  C,. 
determinations  is  known.  Figure  V-l  can  provide  an  estimate 
of  the  number  of  points  needed. 


ERROR 


nnnnrtnnnnnn  'trtn 


Second  -Order  Double-  Interpolation  Scheme. 
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i  0 

2C 

4  0 

3  1 

4* 


DO  10  I=2#NA 
MX  A=  I 
MNAz 1-1 

mma-mxa+i 

IF (HMA.GT.NA  )  MMA=VN/-1 
APxPRBI I .  I  1 
AMODRCMNA.l  I 
AO aPR B ( MM A i 1  I 

IF  I  AP. SE . X. AND. AN.L E .X )  GCTC  25 
COM  IN  'E 
IFLAC=1 
RE  TURN 

00  30  J  -  2  •  NP 
MX  F  •=  j 
MNF= J-l 
MMF  —  MXPH  1 

IF  IMMP.GT.NPI  MMF=MNF-I 
PPsPRfl  f  J  *  ?) 

PN=PRR(MNP,2> 

Pa:ORB(MXP,2) 

IF (PP.GE.T.ANO.PN.LE.Y )  GCTC  45 

CP  NT  IN"<- 

IFLAG=»I 
PF  T'  'BN 

C 1 i =  CP< MNA.MNP > 

Cl  2-CPf  MNA.MXP) 

C  1  3=rCP(  MNA.MMP  ) 

C 2  I -CPI MXA.MNP  ) 

C22-CPI MXA.MXP) 

C  ■>2-FP(MXA,MMP) 

C  I  -  r  PI  MM  A,  MNP  ) 

C 12-CPT MMA.MXP) 
r"  cpi  wma.mmp) 

7  |  -<x-AP>A(X-AO>/(AN-AF»/(AN-AC> 
t  (  x-AN)*IX-AD)/(AP-AN>/(  AP- AC  ) 

r  >;(  *- AN > *(  X -AP ) / < AC- AN )/ < AO- AD J 

1  - r  I • r I  1  ♦  F  2  AC  ?  I »F  34C 1  I 
r  vt  |*ri?AF?*c??FF34C32 
C3T  I  4T  1  3AF24C23  AF34C33 
7  -  ( V  —  P P I • ( V  —  PO I / ( P N — F F  I/IPN-PO  )4CI 
I  { v-PN| • IT-PO ) / (PO-FN l/IPP-PQ  )  ♦  C2 
I y-PN 1 • ( V-PP 1/IPO-f N  )/ I PC-PP »*C3 
tr  ’  irn 
c\r* 


OCCC447C 
OCCCAAec 
^  C  C  44S  C 
r'CfM'OC 
OCCC A5 1 C 

ANOI  PrcrA«-?C 

ANCI  0CCf4f 3C 
OCCC 4540 
C  C  C  C  45*  C 
CCCCASf C 
0CCC457C 
0CCC45F 0 
i np'  ccec45<;e 
KNEOC  f C AECO 
OCCCAE 10 
0CCC4F20 
0CCC4C7C 
OCCCAf 40 
0  C  CC  4E  c  C 
0CC04F 55 
3CCC4C5F 
0CCC4CFC 
CCCC4F7C 
0  C  C  0  4F  7r 
0CCC4FEC 
OCOfAf E C 
TCCC47CC 
OOCC47IO 
0CC0473C 
OCCC474C 
0  C  C  C  4  7*  0 
0CCC475-" 
0CCC47ce 
0CCCA7eC 
0CC0477O 
CCC0477' 
OCCC47EO 
0  C  C  C4  7C  C 
OCCOAPCC 
0CCC4EI0 
00CC4C2C 
CCCC4P4C 
orr C4F45 
OCCCAE'O 
Off  CAFF  C 
0CCC4E70 
CCCC46C0 
occc4e«-c 
0CCC4ECC 
OCCCAF 1 0 
Of CC4F2C 
0CCC4F3T 
0CCC4F4C 
Of CC4F«0 
0"CC4FF0 

♦  CCCCAS7C 

♦  CCCC4FPC 

OCCCARSO 
OCCC5COC 
Cf  CC*0 I f 


SLERCUTINE  C PC AL (NA, HP, FPe. CP. X,  Y, 2.  I  FLAG) 

NA.NP  =  »  OF  ALPHA  ANC  PHI  ANGLES  IN  THE  CP  C AL  I PF AT  I r N 
AINA)  =  VALUES  OF  THE  ALPHAS  IN  THE  CALIBRATION  T API  F  (VAN 
P(NP)  =  VALIJFS  CF  THE  PHIS  IN  THE  CAt.  IBRATICN  T  AEI.  F  (PITCH 
CP(NA.NP)  =  VALUE  CF  CP  F  C  R  EACH  ANCLE  SET  (  AINA).p(NP)  ) 

X  =  DESIRED  ALPHA  ANGLE 

Y  =  DESIRED  PHI  ANGLE 

7  *  CALCULATED  CF 

I  El.  AG  =  ERROR  FL  AC 

THIS  PROGRAM  ESTIMATES  THE  VALUE  CE  CP  FOR  A  GIVEN  ANGULAR 
ALPHA.  PHI.  USING  A  t  INEAR  OCLEL5  I N T EP F CL A T I  ON  SCHEMF  PEUlfFK 
VALUES  OF  CP  FOR  ANGLES  ABC VE  AND  EEIC*  THE  DESIRED  ANGIE. 
DIMENSION  PReiNA.2)  .CPINA.NP) 
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DISTRIBUTION  LIST 


De tense  Ducumen  t  a  t  i  on  Center 
Cdiiie  ron  S  to  t  ion 
Alexandria,  Virginia  22314 

Library 
Code  0212 

Naval  Postgraduate  School 
Monterey,  California  93940 

Office  of  Research  Administration 
Code  0  1 2 A 

Naval  Postgraduate  School 
Monterey,  California  93940 

Chai rman 
Code  67 

Department  of  Aeronautics 
Naval  Postgraduate  School 
Monterey,  California  93940 

Director,  Turbo-Propulsion  Laboratory 
Department  of  Aeronautics 
Naval  Postgraduate  School 
Monterey,  California  93940 

Dr.  H.  J.  Mueller 
Research  Administrator 
Code  3I0A 

Naval  Air  Systems  Command 
Navy  Department 
Washington,  D.  C.  20360 

Mr.  Karl  H.  Guttmann 
Code  330C 

Naval  Air  Systems  Command 
Navy  Department 
Washington,  D.  C.  20360 

Mr.  James  R.  Patton,  Jr. 

Power  Program,  Code  473 
Office  of  Naval  Research 
Arlington,  Virginia  22218 

Commanding  Officer 
Naval  Air  Propulsion  Test  Center 
Attn:  Mr.  Vernon  lubosky 

Trenton,  New  Jersey  08628 


30 


1 


1 


1 


1 


lu.  National  Aeronautic;,  and  Space  Administration  I 

Lewis  Research  Center  (Library) 

2100  Bruokpurk  Road 
Cleveland,  Ohio  ^135 

I  I  .  CAG  L i b  t  a  ry  I 

The  Boeiny  Company 
Seattle,  Washinyton  98l2*f 

I  2  .  L  i  b  ra  ry  I 

General  Electric  Company 
Aircraft  Engine  Technology  Division 
DTO  Mai  I  Drop  H43 
Cincinnati,  Ohio  ^5215 

I  3 •  L i b  rary  I 

Pratt  and  Whitney  Aircraft 
Post  Office  Box  2691 
West  Palm  Beach,  Florida  33^02 

lA.  Library  I 

Pratt  and  Whitney  Aircraft 
East  Hartford,  Connecticut  06108 

15.  Chief,  Fan  and  Compressor  Branch  1 

Mail  S top  5"9 

NASA  Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  Ohio  MI35 

16.  Director,  Whittle  Laboratory  I 

Department  of  Engineering 

Cambridge  University 
ENGLAND 

I  /.  P rof .  D.  Ad ler  1 0 

Technion  Israel  Institute  of  Technology 
Department  of  Mechanical  Engineering 
Hai fa  32000 
I  s  r  ae  I 

lb.  Frof.  F.  A.  E.  Breugelmans  ' 

Institut  von  Karman  de  la  Dynamique  des  Fluides 
72  Chausee  de  Waterloo 
I6A0  Rhode-St.  Genese 
Be  I g i urn 
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L  i  b  ro ry 

Air  Research  Mfg.  Corporation 
Division  at  Garrett  Corporation 
LG2  South  36th  Street 
PhoeniA,  Arizona  85 0 3 

Dr.  T  .  0.  Carta 

United  Technologies  Research  Labs 

LOO  Main  St  ree  t 

Hartford,  Connecticut  06108 

Prof.  Jacques  Cheuvin 
Universite  D 1 Ai x-Ma rse i  1  I e 
1  Rue  Honnorot 
Ma  rse  i  I  I  e  ,  Fr an  ce 

Mr.  James  V.  Davis 
Te  ledyne  CAE 
I  330  Laskey  Road 
Toledo,  Ohio  L36OI 

Mr.  Jean  Fabri 
ONERA 

2 3,  Ave .  de  la  Division  Lee  I  ere 
32  Chati I  Ion 
f  ranee 

Prof.  Or.  Iny  Heinz  E.  Gailus 
Lehrstuhl  und  Institut  fur  St  rah  I  an t i ebe  und 
Turbourbei tsmashinen 
Rhe i n . *Wes t f .  Techn.  Hochschule  Aachen 
Temp le rg raben  55 
5100  Aachen,  Germany 

Professor  J.  P.  Gostelow 

School  of  Mechanical  Engineering 

The  New  South  Wales  Institute  of  Technology 

Aus  t  ra I i a 

DR.  Iny.  Hans-J.  Heineman 

DFVLK-AVA 

Bunsens  tasse  1 C 

3L00  Gottingen,  W.  Germany 

Prof .  Ch .  H i rsch 

Vrije  Universiteit  Brussel 

P le i n I aan  2 

1050  Brussel  Is,  Belgium 


28.  I' iot.  J.  P .  Johnston  I 

S  l ori  t or u  University 

he  pur  ti.cn  t  of  Mechanical  Engineering 
Stanford,  California  94305 

2g.  Plot.  Jack  L.  Kerrebrock,  Chairman  ] 

Aeronautics  and  Astronautics  Department 
31*265  Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts  02139 

30.  Dr.  B.  Lekshmi nurayana  1 

Professor  of  Aerospace  Engineering 

The  Pennsylvania  State  University 

233  Hammond  Building 

University  Park,  Pennsylvania  16802 

31.  Mr.  R.  A.  Langworthy  1 

Army  Aviation  Material  Laboratories 

Department  of  the  Army 
Fort  Eustis,  Virginia  23604 

32.  Dr.  A.  A.  Mikolajczak  1 

Pratt  and  Whitney  Aircraft 

Eng i nee  ring  2H 

East  Hartford,  Connecticut  06108 

33-  Prof.  Dr.  L.  G.  Napoli tano  I 

Di  rector 

Institute  of  Aerodynamics 
University  of  Naples 
Vi  ale  C.  Augusto 
80125  Napoli 
Italy 

34.  Prof.  Erik  Nilsson  I 

I  ns  t  i  tu  t  i  onen  for  Stronini  ngsmask  i  n  tekn  i  k 

Chalmers  Tekniska  Hogskola 
Fack,  402  20  Goteborg  5 
Sweden 

35.  Prof.  Gordon  C.  Oates  I 

Department  of  Aeronautics  and  Astronautics 

University  of  Washington 
Seattle,  Washington  98105 

36.  Prof.  Walter  F.  O'Brian  1 

Mechanical  Engineering  Department 

Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  VA  24061 
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j  /.  i '  r . .  I  .  1)1  .  K  .  C'stwil  i  t  sch 

I  c  v  ■ '  f  '•  l  <v!  u  MO  Lh  b  el  >  U  I  e 

> ,  ;„u  13 

•■•r  -  ,  ■  us  trio 

;<t  .  •  .  ••  .  Paran  jee 

li.jb,  i  ! .  j  i  I  i  >.n  Division 

'...■I.  ..ji  Mcii-nout  i  ca  I  Laboratory 

.  !  t).,  ,  !  1'XJ 

.  ,ul-  ic  -  I  / 

I  rni  i  a 

yj .  K .  l  .  i  cdcotk 

School  o!  Meehan i  ca I  Engineering 
f  r  an  tic) J  Institute  of  Technology 
Era",  -iil  j,  Bedford  MK43  OA  L 
ENGLAND 

‘  hi  .  b i Use  A .  Reese 

Director,  Jet  Propulsion  Center 
chow !  ot  Mechanical  Engineering 
Purdue  University 

or,  I  ..,ftte,  li.Jiona  A 790/ 

s 1 .  i  r.  w .  S  ch  I ach ter 

brown,  6ove r i - S u I zer  Turbomachinery  Ltd 
Oe  p  t  .  1  DE 
t  sche  r  Wyss  P latz 
. h -  802  j  Zurich 
Cwi t  ze  r I  and 

'i  J  .  Prof  .  T .  H  .  Ok  i  i  sh  i 

Professor  of  Mechanical  Engineering 
208  Mechanical  Engineering  Building 
Iowa  State  University 
Arnes  ,  I  owa  500 1  1 

A3-  Dr.  Fernando  Sis  to 

Professor  and  Head  of  Mechanical  Engineering  Department 

Stevens  Institute  of  Technology 

Castle  Point,  Hoboken,  Jew  Jersey  07030 

ui .  Leroy  H.  Smith,  Jr. 

Manager,  Compressor  and  Fan  Technology  Operation 

General  Electric  Company 

Aircraft  Engine  Technology  Division 

S  I  u  m.j  i  i  Drop  H  A3 

Cincinnati,  Ohio  45215 
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l<i> .  L  r .  U  .  1  oLuku  I  t 

h'olc:iMjr,  Department  o!  Aerospace  Engineering 
University  of  Cincinnati 
Cincinnati,  Ohio  A522  I 

A6 .  Mr.  P .  Tran,,, 

Manager,  Research  Labs 
Detroit  Diesel  Allison  Divison 
General  Motors 
P.  0.  Box  89A 

Indianapolis,  Indiana  ^6206 

A7.  Prof.  Dr.  W.  Traupel 

Institut  fur  Thermische  Turbontasch i nen 
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